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^ ■ Abstract 
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1 — i ■ A model for the Pomeron at t = is suggested. It is based on the idea 

of a finite sum of ladder diagrams in QCD. Accordingly, the number of s- 
<^> • channel gluon rungs and correspondingly the powers of logarithms in the 

-»«« ■ forward scattering amplitude depends on the phase space (energy) available, 

O i.e. as energy increases, progressively new prongs with additional gluon rungs 

O 1 in the s-channel open. Explicit expressions for the total cross section involv- 

' ing two and three rungs or, alternatively, three and four prongs (with ln 2 (s) 

and ln 3 (s) as highest terms, respectively) are fitted to the proton-proton and 
' proton-antiproton total cross section data in the accelerator region. Both 

$— i ' QCD calculation and fits to the data indicate fast convergence of the series. 

In the fit, two terms (a constant and a logarithmically rising one) almost sat- 
urate the whole series, the ln 2 (s) term being small and the next one, ln 3 (s), 
negligible. Theoretical predictions for the photon-photon total cross section 
are also given. 
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1 Introduction 



It is widely accepted that the Pomeron in QCD corresponds to an infinite sum 
of gluon ladders with Reggeized gluons on the vertical lines (see Fig. 1), result- 
ing ]T|, |3| in the so-called supercritical behavior a t ~ s Qp(0) , a P (0) > 1, where 
ap(0) is the intercept of the Pomeron trajectory. In that approach, the main contri- 
bution to the inelastic amplitude and to the absorptive part of the elastic amplitude 
in the forward direction arises from the multi-Regge kinematics in the limit s — > oo 
and leading logarithmic approximation. In the next-to-leading logarithmic approxi- 
mation (NLLA), corrections require also the contribution from the quasi-multi-Regge 
kinematics 0. Hence, the subenergies between neighboring s-channel gluons must 
be large enough to be in the Regge domain. At finite total energies, this implies 
that the amplitude is represented by a finite sum of N terms M, where N increases 
like Ins, rather than by the solution of the BFKL integral equation [I], EL pi. The 
interest in the first few terms of the series is related to the fact that the energies 
reached by the present accelerators are not high enough to accommodate a large 
number of s-channel gluons that eventually hadronize and give rise to clusters of 
secondary particles. 

The lowest order diagram is that of two-gluon exchange, first considered by Low 
and Nussinov 0. The next order, involving an s-channel gluon rung was studied e.g. 
in the papers [EJ, [7|. The problem of calculating these diagrams is twofold. The first 
one is connected with the nonperturbative contributions to the scattering amplitude 
in the "soft" region. It may be ignored by "freezing" the running coupling constant 
at some fixed value of the momenta transferred and assuming that the forward 
amplitude can be cast by a smooth interpolation to t — 0. More consistently, one 
introduces a nonperturbative model |§ of the gluon propagator valid also in the 
forward direction. The second problem is more technical: as s —>■ oo the number 
of Feynman diagrams that contribute to the leading order rapidly increases and, in 



the coupling, subleading terms coming both from the neglected diagrams and from 
the calculated ones are present. Although functionally the result is always the sum 
of increasing powers of logarithms, the numerical values of the coefficients entering 
the sum is lost unless all diagrams are calculated. 

Conversely, one can expand the "supercritical" Pomeron ~ s a ^ in powers of 
ln(s). Such and expansion is legitimate within the range of active accelerators, 
i.e. near and below the TeV energy region, where fits to total cross sections by a 
power or logarithms are known |J to be equivalent numerically. Moreover, forward 
scattering data (total cross sections and the ratio of the real to the imaginary part 
of the forward scattering amplitude) do not discriminate even between a single and 
quadratic fit in ln(s) to the data. 

Phenomenologically, more information on the nature of the series can be gained if 
the t dependence is also involved. The well-known (diffractive) dip-bump structure of 
the differential cross section can be roughly imitated by the Glauber series, although 
more refined studies within the dipole Pomeron model (DP) (linear behavior in 
ln(s)) [0] show that the relevant series is not just the Glauber one. A generalization 



of the DP model including higher powers of ln(s) was considered in In a 



recent paper [0 the Pomeron was considered as a finite series of ladder diagrams, 
including one gluon rung besides the Low-Nussinov "Born term" and resulting in a 
constant plus logarithmic term in the total cross section. With a sub-leading Regge 
term added, good fits to pp and pp total as well as differential cross section were 



obtained in |I2]|. There is however a substantial difference between our approach 



and that of Ref. [12| or simple decomposition in powers of ln(s), namely that we 



consider the opening channels (in s) as threshold effects, the relevant prongs being 



2 



separated in rapidity by In s 0) being a parameter related to the average subenergy 
in the ladder. Although such an approach inevitably introduces new parameters, 
we consider it more adequate in the framework of the finite-ladder approach. We 
mention these attempts only for the sake of completeness, although we stick to the 
simplest case of t — 0, where there are hopes to have some connection with the QCD 
calculations. 

In Section 2 we consider a new parametrization for total cross sections based on 
the contribution from a finite series of QCD diagrams with relative weights (coeffi- 
cients) and rapidity gaps to be determined from the data. Each set of the diagrams 
is "active" in "its zone", i.e. the parameters should be fitted in each energy interval 
separately and the relevant solutions should match. The matching procedure will 
be similar to that known for the wave functions in quantum mechanics, i.e. we re- 
quire continuity of the total cross section and of its first derivative. In Section 3 we 
present the result of fits to the pp and pp experimental data. Section 4 is devoted to 
a discussion of the truncated series in QCD and to the calculation of the coefficients 
of the powers of Ins. Finally, in Section 5 we will draw our conclusions. 

2 Description of the model 

The Pomeron contribution to the total cross section is represented in the form 



N 



(') =^0(8 -81) Ob? 1 -8), 



(1) 



i=0 



where 



fi = 2 a v LJ > 

3=0 



(2) 
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So is the prong threshold, 9{x) is the step function and L = ln(s). Here and in the 
following, by s and s respectively, sj (1 GeV 2 ) and s / (1 GeV 2 ) is implied. The main 
assumption in Eq. ([!]) is that the widths of the rapidity gaps ln(so) are the same along 
the ladder. The functions fi are polynomials in L of degree i, corresponding to finite 
gluon ladder diagrams in QCD, where each power of the logarithm collects all the 
relevant diagrams. When s increases and reaches a new threshold, a new prong opens 
adding a new power in L. In the energy region between two neighbouring thresholds, 
the corresponding f iy given in Eq. (|l|), is supposed to represent adequately the total 
cross section. 

In Eq. (HD the sum over N is a finite one, since N is proportional to ln(s), where 
s is the present squared cm. energy. Hence this model is quite different from the 
usual approach where, in the limit s — > oo, the infinite sum of the leading logarithmic 
contributions gives rise to an integral equation for the amplitude. 

To make the idea clearer, we describe the mechanism in the case of three gaps 
(two rungs). To remedy the effect of the first threshold and get a smooth behavior 
at low energies, we have included also a Pomeron daughter, going like ~ 1/s in the 
first two gaps with parameters &o and h\ respectively. Then 

f (s) = a 00 + b /s for s < s , (3) 
/i(s) = a 10 + b 1 /s + a n L for s < s < si , (4) 
f 2 (s) = a 20 + a 2 iL + a 22 L 2 for si < s < si . (5) 

By imposing the requirement of continuity (of the cross section and of its first 
derivative) one constrains the parameters. E.g., from the conditions /i(s ) = /o( s o) 
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and f[(s ) = /o(s ) the relations 

h = On s + b , 

a io — a oo — a n m (so) — a ll 
follow. Furthermore, from ^(sq) = /i(sq) and /^(so) = /i( s o) one S e ^ s 

a 2Q = a 22 ln 2 (so) + a 10 + &i(l + ln(so))/«o > 

02i = an - 2a 22 ln(sjj) - 61/ Sq . 
The same procedure can be repeated for any number of gaps. 

In fitting the model to the data, we rely mainly on pp data that extend to the 
highest (accelerator) energies, to which the Pomeron is particularly sensitive. To 
increase the confidence level, pp data were included in the fit as well. To keep 
the number of the free parameters as small as possible and following the successful 



phenomenological approach of Donnachie and Landshoff [I3[|, a single "effective" 
Reggeon trajectory with intercept a (0) will account for nonleading contributions, 
thus leading to the following form for the total cross section: 

a t (s) = al P) (s) + R(s), (6) 

where (s) is given by Eq. ([!]) and R(s) = as"'- - 1-1 (the parameter a is different 
for pp and pp and is considered as an additional free parameter). 

Ideally, one would let free the width of the gap Sq and consequently the number 
of gluon rungs (highest power of L). Although possible, technically this is very 
difficult. Therefore we considered only the cases of two and three rungs and, for 
each of them, we treated free parameter. 



Notice that the values of the parameters depend on the energy range of the 
fitting procedure. For example, the values of the parameters in / if fitted in "its" 
range, i.e. for s < sq, will get modified in f\ with the higher energy data and 
correspondingly higher order diagrams included. 



3 Fits to the pp and pp data 

As a first attempt, only three rapidity gaps, that correspond to two gluon rungs in 
the ladder were considered. Fits to the pp and pp data were performed from yfs = 4 
GeV up to the highest energy Tevatron data (for pp), including all the results from 
there |jl~4 |. The resulting fit is shown in Fig. 2 with a x 2 /d.o.f. ~ 1.71. The values of 
the fitted parameters are quoted in Table 1. Interestingly, the value of so turns out 
to be very close to 144 GeV 2 , i.e. the value for which the energy range considered 
is covered with equal rapidity gaps uniformly. 

Next, we covered the energy span available in the accelerator region by four 
gaps, resulting in 3 gluon rungs and consequently L 3 as the maximal power. After 
the matching procedure, we are left with ten free parameters: first of all Sq, then 
a 00 ,bo,aii,a22,d32, «33, each determined in its range, while the two a's and a (0) 
are fitted in the whole range of the data. The final value for sq turned out to be 
So — 42.5 GeV 2 resulting in a sequence of energy intervals ending at y/s = 1800 
GeV. It is amusing to notice that a search for the phase space region where the 
production amplitude in the multicluster configuration has a maximum leads, with 
the help of cosmic ray data, to an average "subenergy" < s, >~ 44 GeV 2 jl5| , that 
is very near to the value of sq found in the fit. 
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Fig. 3 shows our fit to the pp and pp total cross section data. The values of 
the fitted parameters are quoted in Table 2. The value of the x 2 /d.o.f. is ~ 1.38, 
much better than in the case of two gluon rungs. It is interesting to observe that 
the coefficients in front of the leading logarithms are related roughly by a factor of 
1/10. The value of the effective Reggeon intercept remains rather low, close to 0.45, 



comparable with the value found in Ref. [16 



4 Explicit iterations of BFKL 

From the theoretical point of view, the phenomenological model of Section 2 cor- 
responds to the explicit evaluation in QCD of gluonic ladders with an increasing 
number of s-channel gluons. This correspondence is far from literal since each term 
of the BFKL series takes into account only the dominant logarithm in the limit 
s — > oo. In the following we give concrete expressions for the forward high energy 
scattering amplitudes for photons and hadrons in the form of an expansion in powers 
of large logarithms in the leading logarithmic approximation. 

We start from known results obtained in paper where an explicit expression 
for the total cross section for hadron-hadron scattering has been obtained. In the 
high energy limit, it is convenient to introduce the Mellin transform of the amplitude 

r°° ( s \ ( s \ ~ U} - 1 Im s A(s, t) 



and its inverse 

Inv 

s 2-ki Js-ioo \m 2 , 

The general expression of A(u,t) in the leading logarithmic approximation has the 



MM) i p^/., A(bjt) 
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form: 

<S> a (k,q)F*(k,q) 



A(u,t) = J d 2 k- 



k 2 (q-k) 2 ' 

where $ a (/c, q) and $ 6 (/c, q) (see next equation) are the impact factors of the colliding 
hadrons a and b, obeying the gauge conditions $ J (0, q) = &(q,q) — (j — a,b). 
The quantity F^(k,q) obeys the BFKL equation: 

d 2 k' A(k, k', q)F»(k>, q) - B(k, k', q)F^k, q) 



uF^k,q) = <S> b (k,q)+ 1 J 

with 



2tt (k - k') 2 

q 2 {k - k') 2 + k 2 (q - k'f + k' 2 (q - kf 



2 



k' 2 (q - k') 

,/ x k 2 (q-k) 2 

B(k,k',q)= ,„ ■ - >v> + 



A;' 2 + (jfe' - A;) 2 (g - A:') 2 + (Jfe - A:') 2 ' 
and 

7 = 3—. 

7T 

The strong coupling a s is assumed to be frozen at a suitable scale set, for example, by 
the external particles. The iteration procedure and the reciprocal Mellin transform 
give (besides we put q — 0): 



lm s A(s,0) r 2 ,$ a (A;,0) 



^ (s) = — ; — = J 



d z k- 



(k 2 y 

where 



P=^'nfi) (T) 



7T \m 2 



and the subsequent iterations begin from $o(k) = $ b (A;,0). In the previous inte- 
gral and everywhere in the following, all the momenta are 2-dimensional Euclidean 
vectors, living in the plane transverse to the one formed by the momenta of the 
colliding particles. 
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For the case of photon-photon scattering one has || (see also the references 
quoted there) 



2 9 f k 2 \ 



where 



777/ 



From the integrals^ 



^9 4m 2 
^ =1 + ^ >1 ' 



7T 



2\2 



T(0 - T(t) + T(t) 



Vt 2 + 4t' 2 



1.443 



we conclude 



1 + 6.4— In 



, , (8) 
where cr is a constant and we used the approximated equality (3 x 1.44)/0.67 = 6.4. 

To obtain the cross section of proton-proton scattering, we use the ansatz of 
Ref. Jl7| for the impact factor of a hadron in terms of its form factor F(q 2 ): 



(q,q)=0. 



Here the 2-dimensional Euclidean vector q is related to the 4-dimensional transferred 
momentum Q by the relation Q 2 = —q 2 < 0. Using the formulae given above, we 
obtain, for a simplified but experimentally acceptable form of proton's form factor, 



$o(k) = ak 2 e 



2-bk 2 



(9) 



1 S. Gevorkyan, private communication 
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where a and h are in GeV 2 . It is convenient to define 



2, _ ®n(k) 



then (for n > 1) 

Jo 1 — x\ J Ji x — 1 V x 



and 



/■oo 

a t (s) = 7i dk 2 Mk 2 )J2Mk 
Jo V 



2^ 



(10) 



The integrations can be performed analytically, due to the simple choice of the 
impact factor in Eq. (|]), and the final result is: 



cr t {s) 

1 

3 



2b 



l + 2(ln2)p + 



7T 



— + 2(ln2) 2 
12 v ; 



P 2 + 



^(ln2)+4(ln2) 3 -^C(3) 



P 3 + 



where p is defined in Eq. (^). 



As stressed above, the coefficients of different powers of ln(s/m 2 ) in Eq. (|TT| ) 
refer to the dominant contribution, at asymptotic energies, for each perturbative 
order. In the fit, instead, the Pomeron contribution is determined only from the 
experimental data at high but finite energies. However, we can obtain a rough esti- 
mate of the importance of the subleading contributions by comparing Eq. (|IID with 
the phenomenological fit of the previous Section, in particular with the amplitude 
fs relative to the last gap of the three rungs case. If we assume a commonly used 



value for the strong coupling, a s ~ 0.5 — 0.7 |17| , we must conclude that subleading 
contributions to the QCD Pomeron are important. This finding may be related to 
the fact that energies reached by the present accelerators are not yet asymptotic. 
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5 Conclusions 



Although high quality fits were not the primary goal of the present study, we still 
may conclude that they are compatible with those existing [TB[], and there is still 
room for further improvement. Our main goal instead was to seek for an adequate 
picture of the Pomeron exchange at t — 0. In our opinion, it is neither an infinite 
sum of gluon ladders as in the BFKL approach |, g g, nor its power expansion. 
In fact, the finite series - call it "threshold approach" - considered in this and our 
previous paper || realizes a non trivial dynamical balance between the total reaction 
energy and the subenergies equally partitioned between the multiperipheral ladders. 

An important finding of the present paper, confirmed both by the fits to the pp 
and pp cross sections (Sec. 3) and by the QCD calculations (Sec. 4) is the rapid 
convergence of the series. 

"Footprints" of the prongs at low energies are slightly visible in Fig. 2 (especially 
in the case of pp scattering where the contribution from secondary Reggeons is 
smaller than in pp). A more detailed study of this phenomenon may answer the 
question whether this is merely an artifact or a manifestation of Pomeron's basic 
properties. The quality of the present fits will be definitely improved when future 
data from RHIC and LHC will be available. As to the QCD calculations, their 
precision and efficiency are biased at finite energy by difficulties in estimating non- 
asymptotic terms, neglected here and elsewhere. As a consequence, e.g. in the 
diagram giving a ln 2 (s) contribution, the constant and ln(s) terms (see Fig. (1)) are 
neglected. 

In any case, few terms of powers of ln(s) are sufficient to describe the data at all 
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realistic energies. As a guess, we do not exclude that higher terms cancel completely, 
but anyway they are negligibly small. The case of two terms (logarithmic rise in 
s) is particularly interesting as it corresponds to a dipole Pomeron with a number 
of attractive features |10| such as self-reproducibility with respect of unitarity cor- 



rections. In case of a ln 2 (s) rise (three terms) we still should not worry about the 
Froissart bound, so ultimately the Pomeron as viewed in this paper does not need 
to be unitarized. This conclusion is an important by-product of our paper. For the 
dipole Pomeron, relevant calculations for t 7^ are interesting and important but 



difficult. In the case of a single gluon rung they were performed in Ref. |T2| and, 
with a non-perturbative gluon propagator, in the last reference of ||. 

The role and the value of the width of the gap, sq, is an important physical 
parameter per se, independent of the model presented above. We have fitted it 
and compared successfully with the prediction from cosmic-ray data. However its 
value may be estimated e.g. as the lowest energy where the Pomeron exchange is 
manifest, although the latter is also a matter of debate. Further fits of the model 
to new experimental data may settle some details left open by this paper. 
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TABLES 



so 


a(0) 


a pp 




147.97(93) 


0.441(11) 


71.7(3.4) 


0.00(37) 


bo 




a u 


0-22 


35.6(1.5) 


38.097(23) 


2.300(38) 


0.857(61) 


a 10 




021 


X 2 /d.ol. 


24.3 


110.1 


-14.85 


1.71 



Table 1: Value of the parameters in the case of 2 rungs. The parameters bi and a... 
are given in units of 1 mb. The quantities in round parenthesis represent the errors. 
Parameters without error are derived from the matching condition. 



s 


a(0) 






bo 


a 00 


42.43(53) 


0.4295(95) 


160.9(8.6) 


85.2(6.6) 


-180(18) 


33.20(29) 


an 


022 


CI32 


^33 


aio 




2.631(86) 


0.324(55) 


0.19(22) 


0.000(38) 


20.7 


38.6 






031 




X 2 /d-oi. 




-2.19 


22.3 


0.72 


-68.42 


1.38 





Table 2: Value of the parameters in the case of 3 rungs. The parameters fej and a 
are given in units of 1 mb. The quantities in round parenthesis represent the errors. 
Parameters without error are derived from the matching condition. 
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Figure 1: Schematic representation of the total cross section in the leading ln(s) 
approximation (first row). Double lines represent protons or anti-protons, vertical 
zig-zag lines are Reggeized gluons, horizontal wavy lines are gluons. The effective 
vertex for two Reggeized gluons and one gluon is defined in the second row. Here 
external lines here can represent quark or gluons. 
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Figure 2: Total cross section calculated up to 2 gluon rungs and fitted to the pp and 
pp data. 
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Figure 3: Total cross section calculated up to 3 gluon rungs and fitted to the pp and 
pp data. 
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